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Abstract 

We show that almost all metric-afhne theories of gravity yield Einstein equations with a 
non-null cosmological constant A. Under certain circumstances and for any dimension, it is also 
possible to incorporate a Weyl vector field and therefore the presence of an anisotropy. The 
viability of these field equations is discussed in view of recent astrophysical observations. 

PACS number: 0450 



Recent observational evidence, pointing to a non-vanishing value of the cosmological constant [1] 
and to the existence of a preferred direction [2] in the Universe, has reopened the debate on which 
is the 'most correct' theory of gravitation. As in any field theory there are two aspects to be 
considered: the field equations and the variational formulation of them. Undoubtedly, general rel- 
ativity is the best theory we have for the description of gravitational phenomena. The Einstein 
field equations were obtained on a phenomenological basis and therefore they exhibit a strong ob- 
servational agreement. On the other hand, the variational formulation of Einstein field equations 
provides a systematic way to deal with them. However, this formulation must not be considered 
as fundamental as the field equations themselves and, in fact, several of the problems of general 
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relativity come from its variational formulation. Just to mention a well known problem, the 'un- 
physical' non-renormalizability of general relativity is obtained from a dimensional analysis of the 
Einstein-Hilbert Lagrangian even when the Einstein field equations describe 'physically' sensible 
situations. Therefore, in order to accommodate the new observational data, not describable in the 
context of general relativity, it would be better to start by looking for field equations with a better 
observational agreement than just on the basis of theoretical considerations. Of course, many ad 
hoc field equations can be conceived to agree with observation and in this case a phenomenological 
argument, such as that leading from Newtonian gravity to general relativity, is lacking. In the face 
of this situation, and in spite of the comments above, we must allow our search for the needed 
modifications to be guided with the help of theoretical considerations. 

It has been shown recently [3] that several families of gravitational Lagrangians give rise to 
the same field equations, namely Einstein field equations with a non-null cosmological constant 
A. Therefore, these field equations are better suited to accommodate the recent observations con- 
cerning a non-null value of A. This 'universality' property was originally observed for Lagrangians 
depending in an arbitrary way from (R) = g^ v ' Rp U (T) [18, 19, 20, 21]. Later on this property was 
established for a wider family of Lagrangians depending at most on (R 2 ) = Rp V {T)R tJ - v '(T) [3]. 

In the present work, we show that these universal field equations can be obtained from a family 
of variational principles which is even wider than that considered previously. In addition, almost 
any metric-affine Lagrangian will serve the prupose of providing a variational description of them. 
Furthermore, we show that, under certain circumstances, it is also possible to incorporate an 
anisotropy through the appearance of a Weyl vector related to conformal invariance. In previous 
works this possibility was established only for n = 2, where n is the dimension of the base space. 
Here we show that this can be done in any dimension. This result is particularly relevant for n = 4 
since a Weyl vector may describe the recently observed anisotropy of the Universe. 

The freedom in the choice of the Lagrangian provides us with a theoretical setting to incorporate 
almost any other desirable property of a viable gravitational theory, in particular renormalizability 
and scale invariance, without restricting the field equations. 

Let us start with a brief reminder of the state of the art in canonical gravity. A large family of 
different situations in gravitational physics is correctly described by the Einstein field equations. 
The problem of obtaining Einstein field equations from a variational principle was first addressed, 
almost contemporarily, by Einstein [4] and by Hilbert [5]; cf. [6] for an interesting historical account. 
Under a suitable set of simplifying assumptions they concluded that a possible Lagrangian for the 
Einstein field equations is given by 

C EH = R( E )g 1 / 2 , (1) 

where -R(g) = g^R^uis) is the Ricci scalar constructed from a purely metric Ricci tensor Rp U (g). 
Almost inmediately, Palatini [7] presented an alternative Lagrangian yielding the same field equa- 
tions, namely 

Cp = g^R^(T)g 1 / 2 , (2) 

where Rp U (Y) is the symmetric part of the Ricci tensor for a connection T x fll/ symmetric in its two 
lower indices. 

For later convenience we will display our convention for the Ricci tensor. The starting point is 
the Riemann tensor given by 

R /o/ii/(r") = OfjY V p d v Y + r r V p — r vc r pp , (3) 
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and its is obviously antisymmetric in its last two indices. From contractions of the Riemann 
tensor alone, not involving any auxiliary field, two second-rank tensors can be obtained, namely 
Ffiu = R X \fj,u = d^u ~ <9^r M , with T M = r A MA , and the Ricci tensor 

iV(r) = R\xu(r) = d x r x ^ - ^r A v + rV r% - r\ a r% . (4) 

In the rest of the work we will restrict our considerations to the symmetric part of the Ricci tensor 
(4) for a connection symmetric in its two lower indices. 

The Lagrangians (1) and (2) both lead to the same set of field equations, namely, the Einstein 
field equations 

iV(g)-7^(g)^ = 0- (5) 

After considering the trace of equation (5), which implies R(g) = 0, the Einstein field equations 
can be rewritten in the simplified, but equivalent form 

iV(g)=0. (6) 

The exception is the dimension n = 2 in which case (5) is an identity. Equations (6) have a strong 
phenomenological support since, in the non-relativistic limit, they are equivalent to Newtonian 
gravity and furthermore they have been tested and found to be correct in a variety of situations 
beyond the realm of Newtonian gravity [8]. However, due to several reasons stemming from the 
need to obtain better observational agreements [9] to formal problems [10], such as those found in 
quantum gravity, alternative theories need to be seriously considered. 

The first modification of generqal relativity is due to Einstein himself [11] and involves the 
introduction of a cosmological constant A. In this case, the Lagrangians (1) and (2) are modified 
to 

C EHA = [R(g)-A n ]g^ 2 , (7) 

and 

£p A = [^V(r)-A n ] 5 1 /2 ) (8) 

where n, as in the rest of this work, denotes the dimension of the base space. In both cases the 
modified field equations are 

R^{g)-\[R{g)-K]g^ = U- (9) 

Considering the trace of this equation, which implies R(g) = nA n /(n — 2), the Einstein field 
equations, with a cosmological constant, can be rewritten in the simplified form 

V(g)- ^-T2) A «V = 0' (10) 

The exception is the case for n = 2 in which the first two terms in (9) are an identity and therefore 
A2 = 0. In view of this result it is more convenient to write A n = (n — 2)A. With this new 
parametrization the field equations (10) are rewritten as 

V(s)- A ^ = o, (11) 
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and the trace is given simply by R(g) = nA. Our later generalizations will be based on this 
simplified, but equivalent, form of Einstein field equations. 

As mentioned above, recent observations [1] point to a non-null value of the cosmological 
constant A and therefore (11), rather than (6), need to be considered as the correct field equations 
for the gravitational field. A small, but non-null, value of the cosmological constant is most welcome 
since it is able to solve several observational problems, such as conflicting ages of stars compared 
with the age of the Universe [12]. The Newtonian limit remains unmodified since A is small enough 
to play a significant role at that regime. 

The Lagrangians (7) and (8) are correct classically, since they lead to the correct field equations. 
However, due to several reasons, mainly non-renormalizability, they do not serve as a starting point 
for a quantization program. As remarked by Martellini [13], what is really relevant for a variational 
formulation of a gravitational theory is to obtain the correct classical equations, in this case (11), 
but the Lagrangian is completely unrestricted for the rest and may even strongly differ from the 
simple prescriptions contained in (7) and (8). This remark is most important when one looks for 
better starting points for a quantum gravity programme. 

The simplest model pointing to satisfy the renormalizability requirements necessary for quantum 
gravity is the Lagrangian R 2 . In the purely metric case we have 

£ RHs) = [aR 2 (g)+PR, u (g)R tlu (g)+7R^x P (g)R tluXp (s)]9 1/2 ■ (12) 

This Lagrangian has been extensively considered from the point of view of necessary counterterms 
for achieving the renormalizability of the Einstein-Hilbert Lagrangian. However, the Lagrangian 
R 2 can be considered by itself, without any reference to general relativity, avoiding in this way 
the introduction of dimensionful coupling constants [14]. In this case it is possible to establish 
an equivalence between (12) and general relativity coupled to some kind of matter. However, the 
Lagrangian (12) is non-polynomial in the metric field and this makes it useless for a quantization 
programme [13]. 

The situation is much better for a metric-affine theory where the starting point is the Lagrangian 

£ RHr) = [ a {R) 2 + (3 (R 2 ) + 7 Rie 2 ] g 1 / 2 ; (13) 

(R) = g^ u R^ u (T), (R 2 ) = Rpv^Wfi), and Rie 2 denotes all possible, algebraically independent, 
contractions of the Riemann tensor R x ppv {T) and since it has fewer symmetries than the Riemann- 
Christoflel tensor R\pnv{g), more terms of that kind need to be considered in (13). In this case 
it is possible to establish the equivalence between (13) and general relativity with a cosmological 
constant, as we will shown below. The Lagrangian (13) is far more appealing than (12) since it is 
polynomial in T and therefore more amenable for a quantum treatment [13]. The terms Rie 2 are 
Yang-Mills-like and they can be dealt with in that context, see [15]. The terms containing (R) and 
(R 2 ) only, were first considered by Stephenson [16] and by Higgs [17]. 

Lagrangians with a more general dependence on (R) and (R 2 ) have appeared in several studies 
related to the gravitational field equations. A generic Lagrangian of the form £^ = f({R))g 1 / 2 

was considered by Lim [18] and, under a simple assumption, he obtained Einstein field equations 
with a cosmological constant, namely (11). This was a first indication that one can consider a 
wider family of variational principles from which to obtain gravitational equations. For a restricted 
family of Lagrangians £(^ R y Hamity and Barraco [19] obtained similar field equations in which, 
in addition to the cosmological constant A, conformal invariance was also incorporated. Finally, 
Ferraris et al. [20, 21] has established the important result that all Lagrangians £^ possess 
'universal' field equations corresponding to Einstein field equations with a cosmological constant A, 
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namely (11). Recently, the universality of the field equations (11) was also proved for Lagrangians 
with an arbitary dependence on (R 2 ) [3]. 

These results could be interpreted as an indication that the cosmological constant and other 
desirable properties of the gravitational theory may be obtained as properties of the field equations, 
rather than being incorporated into the Lagrangian. 

In dimension n = 2, it is possible to incorporate a Weyl vector related to conformal invari- 
ance. The presence of a vector field indicates the existence of a preferred direction. For n = 4 a 
vector field would be a most welcome ingredient in the field equations since this might describe the 
recently observed anisotropy of the universe [2] . 

In the present work we show that the universal field equations (11) can be obtained from a 
wider family of Lagrangians given by C = L((R))g 1 / 2 , with an arbitrary dependence on (R)^ u = 
g^ x Rxu{T). Therefore, it may depend on (R) and (R 2 ), but also on higher-order powers of (R) u ' v . 
More importantly, we show that for any dimension, in particular for n = 4, under certain circum- 
stances, it is possible to incorporate a Weyl vector field W^. This is also possible in the situations 
studied previously by other authors, however, they did not mention this possibility. The Weyl 
vector field W M contributes new terms to the field equations which finally read 

R^u(s) ~\{n-2) (V„W V + V„W„) + \ (n - 2) W„ W v 

-Q(VW) + i(n-2)W 2 ) g^-Ag^ = 0, (14) 

where (VW) = g^V^Wv and W 2 = g^W^Wu- For the special case n = 4 the field equations 
above reduce to 

^V(g) - \ {VpW v + V„W„) + \w„W v -^ ({VW) + W 2 ) 9llv - k 9tiv = . (15) 

Our method of proof is essentially different from the one used in [3, 20, 21] and is strongly based 
on a systematic use of the Cayley-Hamilton theorem, which we will now describe. 

Let us review some fundamental results related to the Cayley-Hamilton theorem. Let us start 
by considering a generic n x n matrix (R)^ ' v . Higher powers of this matrix are denoted by 

(RY „ = (Rfo ■ ■ ■ (R)\ ■ (16) 

k times 

Accordingly, we denote the traces by 

(R k ) = (R)% ■ (17) 

Let us introduce next the cumulant coefficients (C k ) given by 



<C°> = 1, 

(C 1 ) = (R), 

(C 2 ) = \ ((Rf - (R 2 )) , 

(C 3 ) = i (( J R) 3 -3( J R)( J R 2 ) + 2( J R 3 )) , 

(C 4 ) = I; ((R) 4 - 6 (R) 2 (R 2 ) + 8 (R) (R 3 ) + 3 (R 2 ) 2 - 6 (,R 4 )) , 
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(18) 



Essentially the cumulants relate the coefficients in a Taylor series with those of a Fourier expansion, 
some of the cumulants (18) are familiar objects in matrix calculus; namely, (C 1 ) = tr((i?) M ^) and 
(C n ) = det {{R)^ v )- A further property of the cumulants is that, for a given n, only the first n 
cumulants are non-trivial, while upper cumulants are identically zero, (C k ) = 0, k > n. 

The Cayley-Hamilton theorem [22] states that only the first (n— 1) powers of (-R) M „ are linearly 
independent. The relation among them is the characteristic polynomial 

n 
k=0 

Accordingly, only the first n traces (17) are algebraically independent. 

For the purposes of this work we introduce a further set of relevant quantities which are obtained 
by considering derivatives of the cumulants (C h ) as follows. Let 



The first (C k ) , s are given by 



The equation 



{c 2 r v = (R)(Rr u -(Ry v , 

= (c 1 ) (cy v - (cy x (&)\ , 

(CY U = \ ((R) 2 - (R 2 )) {R)\ + ((Ry x - (R) (R)\) (R)\ 

= {c 2 ){cy u -{cy x {cy u , 

'■■ ■ (21) 

(C n r v = (C n )6Z, (22) 

is particularly interesting and is a statement similar to (19) but in terms of cumulants. The product 
of two (C)'s follows the rule 

<cTa (c j )\ = E [(c i+j - (e+1) ) {c e+ y u - (cy x (c i+ j-y u ] , (23) 

where £ = min(i,j). Another interesting relation is given by 

^p. = {C*- 1 ) g»» - (C k ~y v . (24) 
The relations (21) can be inverted to express (-R)'s in terms of (C)'s. Relations (18) are inverted 

to 
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(R) = (C 1 ), 

(R 2 ) = (C 1 f-2(C 2 ), 

(i? 3 ) = (C 1 ) 3 -3(C 1 )(C 2 ) + 3(C 3 ), 

: , (25) 

while relations (21) to 



(rY u = {c i ){c i r v -{cy v , 

(RYu = ((c 1 ) 2 - (c 2 )) {c l T v -{c l ){c 2 r v + {c*T v , 

'■■ • (26) 

Finally we can write the characteristic polynomial (19) in terms of (C)'s just to obtain (22), showing 
that this is, in fact, a statement equivalent to (19) in terms of the cumulants (C)'s. 

In view of the above results we are now ready to write the most general Lagrangian with 
an arbitrary dependence on (R) fl v . We must consider Lagrangians with a dependence restricted 
only to the first n traces {R k ), k = 1, ■ ■ ■ , n and this is equivalently done in terms of the (C fc )'s, 
k = 1, • ■ ■ , n. Namely 

£ n = L n ((C 1 ), ■•■,{C n ))g 1/2 . (27) 

Now we must consider a Palatini-like variational principle in which the metric and the connection 
are varied independently. 

Variation with respect to the metric gives 

(C-Y„ 3$=% + - + (CT„ + (<C> ^ - \ L.) « = 0, (28) 

where we have used equation (22). In order to completely determine the tensor (R)^ u we need to 
know its first (n — 1) powers and for this purpose it is necessary to have (n — 2) other relations 
similar to (28). We may obtain them by contracting equation (28) with {C k )^ v , k = 1,- • • ,n — 2 
and using (22) and (23) to reduce higher powers of C's. We obtain a set of n relations of the form 

/ D\ 
\D n -\ 

where the coefficients -D's and £J's are functions of (C)'s. 

Some exceptional cases may appear. The first one is given by the condition det(-D) = 0, 
but this is satisfied only for a restricted family of Lagrangians with null measure in the space of 
functions. Another exceptional situation is obtained when equation (28) is an identity. In this case, 
no algebraic restrictions are imposed over the Ricci tensor. The solution is 

£e = ^{C^g 112 = ^/det[i^(r)] . (30) 



D l n -i \ ( {Cy v \ ( E 1 



D n - v n -i) \(c n -y u 



U U 1 



(29) 



n—l 
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This is the Eddington Lagrangian [23] which does not depend on the metric, in agreement with the 
fact that (28) is satisfied identically. 

The regular case is characterized by det(-D) / and we arrive at the relations 

(C k r v = (F k )5», (31) 

with (-F)'s functions of (C)'s. Therefore, under quite general assumptions the field equations (28) 
are equivalent to (31). However, this result can be simplified even more since for k = 1, equation 
(31) reduces to 

(Cy^iRy^iF^SZ. (32) 

The cumulants (18) then reduce to 

/r k \= nl (f 1 )* = nl 

K ' (n-k)\k\ K ' k\{n-k)\ 
Therefore, equation (32) is rewritten as 

R uv (T) = F{(R))g llv , (34) 

where F {(R)) is some function of (R). The trace of this equation gives (R) = nF((R)). Two 
important situations may appear: 

1. If (R)= constant, then equation (34) reduces to 



n 



(R) 



(33) 



iV(r)-A«7^ = 0, (35) 

with A = constant. Equation (35) above looks the same as the Einstein field equations (11). 
However, they are not yet equivalent since a metricity condition, i.e., a relationship between 
r and g, is still required. This relationship will be obtained starting from equation (46). 

2. The second situation appears when the trace of equation (34) is identically zero. This means 
that (34) must be of the form 

iV(r)-^(i?)v = °- ( 36 ) 

This last situation is the most important for physical applications and we will return to it 
below. 

As a way to verify the correctness of the procedure developed above, let us consider in detail 
the case n = 3. The Lagrangian is given by 

C 3 = L 3 ((C 1 ), (C 2 ), (C 3 )) g 1 ' 2 . (37) 
Variation with respect to the metric gives 

< cl )".^ + ^ + (< c3 > wk - H * - o. < 38 > 

where we have used (22) for n = 3. The trace of this equation gives 
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<C 3 > 



dL 3 1 

a(c 3 ) 2 3 



i 

3 



2(C 2 ) 



dL 3 

d(c 2 } 



+ {&) 



dL 3 



Therefore, equation (38) is rewritten as 



i 

3 



2(C 2 ) 



d{c 2 ) 



2 dL 3 + dL 3 



d{c l ) ' 1 v d(c 2 ) 

Let us now contract this equation with (C 1 )^ v . Using (23) we then obtain 

dL 



6Z = 0. 



{ C 2 )^- + 2(C 1 ) 



d(c 2 ) 



die 1 ) 



^ - {c 2 r„- (c 3 ) =o. 



d(c 2 ) v 



Equations (40) and (41) allow to solve for (C 1 )^ and (C 2 )^, as in (31), if 



1 dU 



3 d(C 2 } 



(C 2 



dL 3 

d(c 2 ) 



+ 2(C 1 ) 



die 1 ) 



+ 



( dL 3 



The condition A = is fulfilled only by a restricted family of Lagrangians, namely, 



L 3 = <CV/3(4„- l)Ve 



±1/3 ( i + 2vr^ y 1/6 



(39) 



(40) 



(41) 



(42) 



(43) 



where u = (C 2 ) / (C 1 ) 2 . Only one sign can be chosen in (43) and G is determined from equation 
(39). This shows that our derivation is correct, except for a null-measure set. 

Let us now turn to the second set of field equations, the metricity condition, which must be 
obtained from (27). Variation with respect to the connection V gives 



V A ( 



7"V /2 )=0, 



(44) 



with 



y 1 " = 

dRf+u 



dL r , 
I 



(45) 

where we have used (24). This is the Eddington prescription [23] to define a metric in a metric- 
affine, or purely affine, theory. Let us now consider the two different situations which appeared 
from (34). 

1. Combining equation (45) with equation (35) we obtain 7^ = Kg^ u , where K = constant. 
Therefore, g^ satisfies equation (44), with 7^ = g^ v . The solution to (44) is 

r V = { A } (s) = \ 9 Xp ( d »9u P + d v g w - d p9tlu ) (46) 

and corresponds to the Christoffel symbol for the metric tensor g^ u . Replacing (46) into (35) 
we obtain the universal field equations (11). The exception is n = 2 and this case can be 
considered under the following point. 
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2. The second case is much more interesting. In this case the relevant equation is (36). This 
equation implies that there is a fundamental invariance in the theory. In fact, nothing changes 
under the transformation g^ u — ► Q^xjg^. This corresponds to conformal invariance and all 
of the relevant equations should display this invariance property. Let us remind ourselves 
that g 1 / 2 — > iW 2 ^ 1 / 2 . Therefore, conformal invariance is guaranteed if L n — > Qr n l 2 L n . 
Then 7"" -► ^-^V- The relation (44) can be rewritten in a way displaying this property 
explicitly as 

Va7^"^ (7.rV A f T )7'"' = 0. (47) 
With all the ingredients above we can write ^ as 

7^=(^§) W2 V- (48) 

This construction procedure forces us to introduce, for dimensional reasons, the constant A 
with dimensions dim[A] = L~ 2 ; the factor n is there just for computational convenience. 

As a consequence of conformal invariance, equation (47) is traceless with respect to 7^ and, 
therefore, there is an arbitrariness in the solution of (47). The most general solution of (47), 
incorporating this arbitrariness, is 

Va7^ = -Wa7*u>- (49) 
This equation can be solved for the connection to obtain 

rV = {i}(7) + ^V(7> W), (50) 
where j^j (7) is the Christoffel symbol (46) for the tensor 7^ and 

w V(7, W) = \ ^ (V; lvp + w v 7w - w; 7^) , (si) 

where is an arbitrary Weyl vector field [24]. 

Let us now return to equation (48). This relation can be rewritten as 

lnv = g^u , (52) 
where is a dimensionless function given by 

e * - (m) ■ (53) 

In terms of this parametrization, the connection is now given by 

r A p = {i}(g)+^V(g.w), (54) 

where W l _ l = W' fl + d^. 

Let us finally note that the resulting theory is conformally invariant. One way of breaking the 
conformal invariance is to let ^ acquire a specific value. Due to the conformal invariance, we can 
fix this value at will, in particular, we can choose ^ = 0. In this last case we obtain an effective 
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theory for low energies where conformal invariance has been broken. In this case the field equations 
(36) and (53), evaluated for the connection (54) with ^ = are 



(n-2) 



(V M W„ + V„Wp) + 



(n-2) 



w u w t 



4 



4 




(n-2) 



4 



VF 2 ^ = 0, 



(55) 



R(g) - (n - 1) (W) + 



(n-l)(n-2) 



W 2 = nA, 



(56) 



4 



where now V is the usual covariant derivative with respect to the metric g^ v . The field equations 
(55) and (56) are written completely in terms of the Riemannian metric g^ v and therefore can be 
compared, interpreted and tested in familiar physical situations. 

We can now turn to the conclusions. Firstly, we have shown that almost all Lagrangians depend- 
ing on (R} >1 ' v = g^R\ v (Y), in a generic and almost arbitrary way, possess the same 'universal' field 
equations, namely (11). In the circumstances described under point (2), it is furthermore possible 
to incorporate a preferred direction. The freedom contained in the choice of the starting Lagrangian 
may be used to satisfy the several different requirements usually needed in order to reach a consis- 
tent quantum version of canonical gravity. For n = 4 the universal field equations are (15) and they 
correspond to the usual Einstein field equations with a necessarily non-null cosmological constant 
and a Weyl vector field W M . Both of these new properties of the gravitational field equations have 
strong observational support. Naturally, the next step is to look for solutions of equation (15) 
in particularly interesting physical situations. A and Wu are of cosmological character and their 
values must be determined from astrophysical observations. Let us first analyse the case = 0. 
Recent observations [1] point to a non-null (even when small) value of the cosmological constant A. 
A non-null cosmological constant does not have a significant role over solar system scales. In fact, 
the corresponding spherically symmetric solution is the Kottler metric. The Schwarzschild solution 
is naturally recovered for A = and for small A, which is the actual case, there are no observational 
effects. However, observational effects may be expected at the level of galactic dynamics [25]. On 
the other hand, there is some evidence that the universe is not completely isotropic, as assumed in 
the cosmological principle, and recent observations point to the existence of a preferred direction in 
the Universe [2]. That direction may be accounted for by the Weyl vector field appearing in (15). 
The possibility of accommodating the observed anisotropy in the field equations (15) is currently 
under study with the use of Bianchi spaces (the natural anisotropic extension of FRW spaces). 
Once the value of lf„ is determined, we can look for other observational effects of W^. However, 
and in a way similar to the case of the cosmological constant, we believe that the magnitude of 
the anisotropy is small as to play a significant role over solar system scales. Therefore, observable 
effects must be sought at the level of galactic dynamics. 

In conclusion, the 'universal field equations' (15) can be considered as a viable generalization of 
the Einstein field equations of general relativity. Furthermore, almost any metric-afnne Lagrangian 
is a good candidate for this purpose and this additional freedom may serve to accommodate several 
other requirements. 

Let us finally add that the universal field equations (15) are not exclusive to the Lagrangians 
considered here. They have also been obtained in [26], in the context of the recently developed 
fourth-rank gravity [27]. 
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